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We study the transition from laminar to chaotic behavior in deterministic chaotic coupled map 
lattices and in an extension of the stochastic Domany-Kinzel cellular automaton For the 

deterministic coupled map lattices we find evidence that "solitons" can change the nature of the 
transition: for short soliton lifetimes it is of second order, while for longer but finite lifetimes, it 
is more reminiscent of a first order transition. In the second order regime the deterministic model 
behaves like Directed Percolation with infinitely many absorbing states; we present evidence obtained 
from the study of bulk properties and the spreading of chaotic seeds in a laminar background. To 
study the influence of the solitons more speciflcally, we introduce a soliton including variant of 
the stochastic Domany-Kinzel cellular automaton. Similar to the deterministic model, we flnd a 
transition from second to first order behavior due to the solitons, both in a mean field analysis and 
in a numerical study of the statistical properties of this stochastic model. Our study illustrates that 
under the appropriate mapping some deterministic chaotic systems behave like stochastic models; 
but it is hard to know precisely which degrees of freedom need to be included in such description. 

PACS numbers: 05.45.+b, 05.70.Jk, 47.27.Cn 



I. INTRODUCTION. 

Spatiotemporal chaos (STC) occurs in many spatially 
extended deterministic systems and remains notoriously 
difficult to characterize Therefore one may attempt 
to map such deterministic chaotic systems onto stochas- 
tic models for which many more analytical methods are 
available. It is then tacitly assumed that, after sufficient 
coarse graining of the deterministic model, the role of de- 
terministic chaos can be taken over by the noise in the 
stochastic system. A critical test of the validity of such 
mappings are the predictions for the transitions between 
qualitatively different states that extended chaotic sys- 
tems display. The key question is then: are transitions 
in deterministic chaotic systems governed by the univer- 
sality classes of stochastic systems? 

As is known for a variety of spatiotemporal chaotic sys- 
tems ^ and as we will show below for the deterministic 
system at hand, chaotic states in extended systems often 
display a mixture of rather regularly propagating struc- 
tures and more disordered behavior. When the propa- 
gating structures, that we will refer to as "solitons" (fol- 
lowing Ql) have a finite lifetime, it may seem that they 
can be ignored after sufficient coarse graining. We will 
find strong indications that this is not always the case, 
and we will give an example where their influence may 
even be so strong as to change the nature of the transi- 
tion. We will also show that extending simple stochastic 
models with the appropriate solitonic degrees of freedom 
can mimic this behavior quite accurately: not only can 
we change the order of the transition, we can also get 
transient non-universal scaling of the type observed in 



coupled map lattices Therefore we conclude that, in 
many cases, deterministic chaotic systems can indeed be 
mapped to stochastic models. A short account of our 
work has already been published |6| . 



A. Historical background 

Chate and Manneville introduced the notion of a 
universal transition to extended chaos via "spatiotempo- 
ral intermittency" (STI) |, | in a study of the deter- 
ministic damped Kuramoto-Sivashinsky partial differen- 
tial equation STI states are composed of "turbu- 
lent" (chaotic) and "laminar" (ordered) patches, and the 
laminar patches remain so except for contamination by 
turbulence at their boundaries. These states are con- 
jectured to occur quite generally when, locally, laminar 
and turbulent dynamics are separated by a subcritical 
bifurcation, and indeed a large number of different ex- 
perimental systems and theoretical models display STI 

As a function of their parameters, STI systems display 
a transition from states where the turbulence eventually 
dies out to states where the turbulence spreads and dom- 
inates. Pomeau proposed jlll an analogy between this 
transition and the phase transition of the stochastic pro- 
cess known as Directed Percolation (DP); for an intro- 
duction to DP, see e.g. jl^, |l^. In directed percolation 
one considers the spreading of "activity" in an absorb- 
ing, inactive background. Earlier, Grassberger and 
Janssen ||T^ had conjectured that any stochastic process 
with an unique absorbing state should be in the same di- 
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rected percolation universality class. 

Relating laminar to inactive and turbulent to active 
states appears to map spatiotemporal intermittency to 
directed percolation. To verify whether deterministic 
chaotic models with an absorbing state would be in the 
DP universality class, Chate and Manneville introduced 
a very simple coupled map lattice (CML) that displays 
STI and numerically obtained the critical exponents that 
characterize the transition from inactive to active states. 
Surprisingly, these critical exponents appear to vary with 
the parameters and are in general different from the DP 
values. Therefore the Chate-Manneville model appears 
to be not in the DP universality class and not even uni- 
versal. 

Grassberger and Schreiber [Q pointed out that the 
presence of long lived traveling structures which they call 
"solitons" in the Chate-Manneville model may lead to 
large crossover times, and conjectured that in the long- 
time limit the behavior of the Chate-Manneville model 
would be in the DP universality class. 

Recently, the Chate-Manneville model with an asyn- 
chronous update rule was studied | [T^ . Here random 
sites are chosen to be iterated forward while keeping the 
others unaltered. For this model, the solitons observed 
for the standard synchronous update rule are suppressed 
and the critical exponents are universal with DP values, 
implying that the Chate-Manneville model with asyn- 
chronous updating belongs to the DP universality class. 
However, the asynchronous updating introduces an ele- 
ment of stochasticity into the model, thus ruining the 
deterministic character of the original model. 



B. Outline 



difficult task to faithfully map a deterministic system to a 
stochastic counterpart. In this particular case, localized 
propagating structures can be identified as responsible for 
the breakdown of DP universality, but one can imagine 
that less easily identified properties of the deterministic 
dynamics could be responsible for such a breakdown in 
other systems. 

The outline of this paper is as follows. In section II we 
discuss the coupled map lattices. Starting from a brief 
discussion of the classic Chate-Manneville model, we in- 
troduce o ur ex tension to lattices of two-dimensional maps 
in section [I A . We show that the new parameter that is 
introduced has a profound effect on the importance of 
"solitons" , and that long living solitons change the tran- 
sition from inactive to active st ates from a second to a 
first order transition in section II B. In the second or- 



der regime, we estimate the bulk critical expo nents using 



finite size scaling techniques in sect ion II C , and mea- 
sure spreading exponents in section [ID , All this data 
is consistent with the coupled map lattice being in the 
universality class of Directed Percolation with infinitely 
many absorbing states, provided that soliton lifetimes are 
short. In section III we discuss the extension of the stan- 



dard Domany-Kinzel cellular automata which includes 
new degrees of freedom that mimic the solitons of the 
coupled map lattices. The mean field equations for this 
model are studied in section HI B, and these show a tran- 



sition from second to first order behavior as a function 
of the soliton lifetimes. We study the phenomenology 
and its statistical bulk proper ties o f the full model in the 
soliton rich regime in section HI C . The behavior of the 
model in the soliton rich regime is quite distinct from an 
ordinary second order transition. 



In this paper we will study a deterministic extension 
of the Chate-Manneville CML that facilitates the tuning 
of the soliton properties. We will demonstrate that the 
influence of solitons may be much more profound than 
setting a crossover time, since they appear to be able to 
change the type of transition from second to first order. 
The role of the solitons is further illustrated in an exten- 
sion of the stochastic Domany-Kinzel cellular automaton. 
In its standard form, all sites of this model can be either 
active or inactive, but we will add a "solitonic" degree of 
freedom that mimics the behavior of the solitons in the 
CML's. The mean field equations of this stochastic model 
show a transition from second order DP-like, behavior to 
a first order transition when the soliton lifetimes are in- 
creased. Numerical studies of this stochastic model also 
find evidence for such a crossover to first order behav- 
ior, although it is very difficult to asses the asymptotic 
behavior for our model. In any case, we present strong 
numerical evidence that the transition is not an ordinary 
second order transition and that there is no asymptotic 
scaling regime, although there are appears to be a tran- 
sient which displays non-universal scaling behavior. 

Our study illustrates that for extended systems it is a 



II. COUPLED MAP LATTICES. 

The model introduced by Chate and Manneville con- 
sists of coupled maps, each of which either performs 
"laminar" or chaotic motion. The model was moti- 
vated by the fact that studies of the deterministic partial 
differential equations, such as the damped Kuramoto- 
Sivashinsky equation, are numerically quite demanding 
and had not provided enough precision to allow a defini- 
tive comparison to DP (^,0]. In one spatial dimension 
their coupled map lattice was defined according to 

u^in + 1) = f{u^{n)) + ^AfU,{n) (1) 

where the subscripts i denote the spatial position, n is the 
discrete time and AfUi{n) ~ f{ui-i{n)) — 2f{ui{n)) + 
/(ui-i-i(n)). This expression is a discrete approximation 
of diffusive coupling in one dimension and introduces spa- 
tial correlations in the system; the parameter e is a mea- 
sure of the coupling strength between a site i and its two 
nearest neighbors at sites (« — 1) and (i -I- 1). 

The map / is chosen such that locally the scalar field Ui 
can be in either of two states: the absorbing (laminar) or 



the chaotic (turbulent) one. When u < 1, / is a standard 



tent map of the form f(u) = 



i|) which dis- 



plays chaotic behavior, while in the region where u > 1, 
/ is simply the identity and leads to a laminar state. 
The sharp discontinuity in / ensures that the two states 
are distinguishable at each site. The parameter r > 2 
determines the steepness of the tent map as well as the 
transition ratio from the chaotic to the laminar regime 
in the absence of coupling. 

The form of the diffusive coupling ensures that turbu- 
lent sites cannot be spontaneously generated in a back- 
ground of laminar sites: states where all sites are laminar 
remain so, and the laminar state is truly absorbing. The 
laminar state is not unique: Updating a state where all 
sites are in the laminar regime {ui > 1) leads, via the dif- 
fusion operator, to a state where all variables are equal 
to the global average value u. 

Once initiated, turbulent activity can spread through 
this CML by infecting laminar patches from their bound- 
aries. The effectiveness of the resulting spreading of the 
chaos depends on the values of r and e. Suppose we study 
the behavior of this system by keeping r fixed while vary- 
ing the coupling strength e. Completely analogous to 
DP, a critical value e — ec{r) exists, such that for e < £c 
an absorbing state is reached with unit probability while 
sustained chaotic behavior (in the thermodynamic limit) 
is found for e > e^- Taking the density of chaotic sites 
or "activity" m as an order parameter, transitions from 
a "laminar" state (where m decays to zero) to a "turbu- 
lent" state (where m reaches a finite value in an infinite 
system) can be studied. 



A. Extensions to two-dimensional maps 

Coupled map lattice's can, in principle, be related to 
continuous time physical systems of weakly coupled ele- 
ments by interpreting the map / as a return map on a 
Poincare section. The time spent by two different sites 
between successive returns would in general be different 
for systems without periodic external forcing, and this 
was precisely the motivation for the asynchronous update 
rule in p^ . However, here we wish to mimic the varia- 
tions in return times in a deterministic fashion. This mo- 
tivated us to introduce a second field in the CML. Note 
that the simplest chaotic oscillator would be a system of 
three phase space dimensions, like the Lorenz equations. 
Applying a Poincare section reduces such a system to a 
two-dimensional map. This is also the case in systems 
with external periodic forcing. Here the simplest real- 
ization would be systems like a damped nonlinear pen- 
dulum or Duffing oscillator with a time-periodic forcing. 
A Poincare section again reduces the system to a two- 
dimensional map and after the synchronous iteration the 
respective units are still at equal time. 

We therefore replace the single variable map f{u) used 
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FIG. 1: Average profiles of a right (a) and left (b) moving soli- 
ton that occurs in the Chate-Manneville model at criticality 
for r = 3. 



in Eq. (^, by a new map with an additional variable v: 

Ui{n+1) = f{ui{n) + -AfU,{n)+Vi{n), (2) 
Vi{n + 1) = b{u^{n + 1) - Ui{n)) . (3) 



Here / is the same map as before and the new parameter 
b is the Jacobian of the full two-dimensional local map; 
this map is invertible for any non-zero b and becomes in- 
creasingly two-dimensional with |&|. The change in the 
local map (^ is analogous to how the 2D Henon map 
|l7| is constructed from the ID logistic map, except that 
b{ui{n +1) ~ Ui{n)) appears here on the right hand side 
instead of bui{n). This ensures that the absorbing state 
fixed points Ui{n) — u* of the old CML (|l|) are mapped 
to the laminar fixed point {ui(n),Vi{n)) = (m*,0) in the 
new CML. The model (|^||) is a completely deterministic 
system with no element of stochasticity and is updated 
synchronously. The value of Ui determines, as in model 
(n^) , whether a given site is " active" or " inactive" . Start- 
ing from the Chate-Manneville case [b — 0) we can follow 
the transition between laminar and chaotic states. As we 
will see below, the new parameter b actually opens up the 
possibility to study the effect of the "solitons" on the dy- 
namical states and transitions of CMLs; this appears to 
be a more important issue than the dimensionality of the 
local map. 



B. Qualitative properties. 

Our CML now contains three freely adjustable param- 
eters (r, £ and 6), and clearly we will have to focus on a 
subset of parameters. Our main focus will be on the case 
where r = 3, although we will also study the transition 
for r — 2.2. For r — i and 6 = the dynamics shows 
many solitons (see Fig. |l|) and the critical exponents ap- 
pear to differ significantly from those of DP. 

To get a feeling for the location of the transition as 
function of b and £, we show in figure |^ the activity (de- 
fined as the average number of active sites) after 1000 
iterations in the ranges —0.3 <b< 0.3 and < e < 0.4 



4 




-0.3 -0.1 0.1 b 0.3 

FIG. 2: Activity in the model (||)-(|) at t = 1000. White 
regions correspond to points in the (e, fe)-plane where the ini- 
tial activity has decayed into an absorbing configuration and 
the darker regions to points with a non-vanishing order pa- 
rameter. Clearly the transition curve becomes quite compli- 
cated; the two branches discussed in this paper are indicated 
as "Branch A" and "Branch B" (see text). The dashed line 
indicates the Chate-Manneville model (6 = 0) 




-2.10"-^ 58 4.10"^2.10"^ 58 

FIG. 3: Activity in the model for 128 systems of size 

L — 2048 and r — 3. (a) Activity as a function of 5e (distance 
to the critical point) at time 2x 10^ for b = —0.1 (squares), b = 
(open circles) and b = 0.2 (closed circles). The transition 
appears much sharper for negative values of b. (b) Steepening 
for the transition at 6 = —0.1 for increasing times: 5 x 10"^ 
(-I-), 5 X 10* (*), 5 X 10^ (x). To stress the magnified scale of 
5e, also the data shown in panel (a) for fe = is plotted (open 
circles) . 



at r — 3.0. The "traditional" transition is that occurring 
at 6 = and e — 0.35984 .... Clearly, for negative val- 
ues of b two additional transitions emerge. Here we only 
study points on the two transition branches labeled "A" 
and "B" in Fig. H; below we focus on the behavior along 
branch A. 
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FIG. 4: Spacetime plots of our coupled map lattice ^ -(g) for 
r — 3 above (left column) and below (right column) criticality. 
Inactive sites are white, chaotic sites are black, (a) b — —0.1, 
£ = 0.353 (b) b = -0.1, e = 0.343 (c) 6 = e = 0.361 (d) 
6 = 0, e = 0.351 (e) b = 0.2, e = 0.374 (f) b = 0.2, e = 0.364. 



1. Qualitative changes in behavior along branch A 

Figure |2| hints that the sharpness of the transition 
varies along branch A: the jump in order parameter ap- 
pears to become steeper for negative values of b. The 
differences in the nature of the transitions are illustrated 
more clearly in Fig. ^ by plotting the value of the order 
parameter as a function of e for 5 = and b = —0.1 av- 
eraged over an ensemble of 32 systems for a number of 
times. The behavior for 6 = is consistent with a contin- 
uous transition, whereas for b = —0.1 longer times lead 
to a marked steepening, consistent with the emergence 
of a discontinuity. 

a. Soliton regime Some effects of the parameter b 
on the dynamics can also be seen from the evolution of 
the binary patterns at r = 3 (Fig. ^). For b = 0, solitons 
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FIG. 5: Examples of good rescaling plots for b = 0.2 
(a,b) and poor rescaling for b = —0.1 (c,d). (a) Absorp- 
tion time r vs. systemsize L, for r — 3,b — 0.2 and 
£ = 0.3727,0.37322,0.37323 (critical value), 0.37325,0.3733 
and 0.3735. (b) Rescaled average activity m' ~ mL~'^ versus 
rescaled time t' ■- t/L^ for r = 3, 6 = 0.2 and e = 0.37323 for 
L — 32, 64, 128, 256 and 512, showing a good data collapse, 
(c) Absorption time r vs. systemsize L, for r — 3,b — —0.1 
and £ = 0.35200, 0.35203, 0.35205, 0.35206 and 0.35207. Even 
small changes in e lead to substantial changes in the absorp- 
tion time, and it is difficult to estimate the critical value of 
e. (d) Rescaled average activity m' versus rescaled time t' for 
r = 3, 6 = -0.1 and e = 0.35203 for L = 64, 128, 256, 512 and 
1024, showing poor data collapse; either the initial decay or 
the tails do not overlap; shown here is a compromise. Note 
that the initial decay is very slow, leading to a small estimate 
for the value of 9. 

can be seen both above and below threshold (Fig. ^-d). 
They consist of pairs of active sites and propagate with 
velocity one. Their maximal lifetimes are of order 100 
(Fig. ^). When b is decreased to a value of —0.1, the 
typical lifetimes of solitons become so long that they typ- 
ically only vanish when they collide with other solitons or 
propagate into turbulent structures. When two solitons 
collide, they either annihilate or create new turbulent 
structures. Such creation is clearly visible in Fig. ^ for 
n « 200 and i w 600. 

For sufficiently large b the isolated solitons present in 
the original model (6 = 0) are suppressed: solitons with 
a lifetime longer than a few iterations are rare here. On 
the other hand, there are regular "edge" states visible, 
where an active state propagates ballistically while emit- 
ting new activity; one example is visible in Fig. ^ for 
n « 400 and i « 800. These structures do not seem to 
influence the order of the transition, but they may very 
well lead to rather large crossover scales. 

In conclusion: the value of b has a large influence on 
the presence of solitons, and also influences the steep- 
ness of the transition. In fact; discontinuities are found 
at points in (e, 6)-space where solitons dominate the dy- 
namics. This implies that the (colliding) solitons have 
a strong influence on the global dynamics and are able 



to change the nature of the transition from a continuous 
to what appears as a first order one. We will make this 
point more precise below. 

C. Finite size scaling in 2nd order regime 

Stochastic systems belonging to the DP universality 
class are characterized by a set of critical exponents de- 
scribing, e.g., the order parameter m(e, L, t) and the be- 
havior of the "absorption time" T(r, e, L), i. e. the av- 
eraged time it takes the system, starting from a random 
initial state, to reach the absorbing state. From finite-size 
scaling arguments [^Sj, one finds that the order param- 
eter m at the critical point Ec should behave as 

m{L,t)^L-f'/''-g{t/L^) . (4) 

For a finite lattice, the absorption time t then increases 
as 

T^L^. (5) 

Finally, for short times (t < L^), g{t/L^) - {t/ U^yP/^w , 
so that for short times m should decay as 

m{L,t) - for t « (6) 

Here the usual dynamical exponent z = i'\\/i'± has been 
introduced, defined as the ratio between the correlation 
length exponent in the time direction and the corre- 
lation length exponent in space v^. The scaling relation 
6 — — /?/^|| connects the critical exponents. 

To estimate the critical exponents for our CML we per- 
formed direct numerical simulations and calculated the 
absorption time r and the order parameter m, defined as 
the average activity. We used ensembles of initial condi- 
tions, in which all initial u values are assigned a random 
value in the chaotic phase, < Mi(0) < 1. The w-values 
of the initial state are set to zero to ensure that they do 
not influence the w-values from the onset of iteration and 
that the analogy with the original model and our variant 
at 6 = is satisfied. 

The behavior of the absorption time at criticality is 
used to determine the critical point and the z exponent 
| ]l6| , |l8|. An ensemble of 128 systems is iterated forward 
in time until an absorbing configuration is reached. The 
average number of time steps needed before reaching such 
a configuration yields the absorption time t. Examples 
of T as function of L are shown in Fig. ||a,c; the best fit 
to a straight slope determines the critical exponent z. 

In Fig. pp,d we plot examples of m' := mL~^ as func- 
tion oit':—t/L^ for a range of L's. When proper scaling 
occurs, as is the case in Fig. |b, the curves for different L 
fall on top of each-other, and the initial power-law decay 
of m' determines the exponent 9. Here an ensemble of 
1000 systems was used. The order parameter was cal- 
culated as the sum of active sites divided by the total 
amount of sites. The systems are iterated forward until 
t « L^, where the algebraic behavior clearly ends. 
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FIG. 6: Spreading of a single turbulent seed through the "natural-initial-state" below criticality, at criticality and above criticality, 
for b — 0.2 (a-c), the Chate-Manneville model for 6 = (d-f) and b — —0.1 (g-i). In all cases r — 3.0. 



Estimates of critical exponents have been done for r = 
2.2 and r = 3.0. For r = 3.0, the critical exponents for 
the original model (6 ~ 0) shows significant deviations 
from the corresponding DP values and the computational 
costs are tolerable. In Fig. |5^,b we show examples of the 
rescaling plots for r — 3,b — 0.2, where a nice data- 
collapse occurs and the transition appears to be of 2nd 
order, and for r = 3, 6 = —0.1 (Fig. ^,d), where the 
data-collapse is poor and the transition appears to be no 
longer continuous. 



The values of the critical exponents are given in Table 
H and correspond simply to the best possible values, irre- 
spective of the quality of the data collapse. For r ~ 2.2 
DP values are found for |6| > 0.01. For r = 3, the crit- 
ical transition on Branch B appear to be DP-like, while 
on Branch A a crossover to DP values is found when b 
is large enough ( \b\ > 0.15). This regime coincides with 
values of b where the solitons are suppressed in the space- 
time plots, and a continuous transition takes place. The 
soliton dominated dynamics at & = —0.1 is reflected in 
the extremely low value of the exponent 9, characterizing 
the decay of the order parameter. Here the data-collapse 
is rather poor as shown in Fig. ^,d. 



D. Spreading of Turbulence in 2nd order regime. 

So far the critical properties of the CML's starting from 
"homogeneous" states have been studied, i.e., with ini- 
tial conditions where each site in the lattice is assigned 
a random number in the chaotic (turbulent) phase. A 
different approach is to consider the spreading of a single 
turbulent seed in an otherwise laminar configuration (see 
Fig. ^). This makes it possible to study the dynamical 
critical exponents, or spreading exponents, and see how 
these compare to the directed percolation counterparts. 

For spreading of activity in stochastic systems with 
absorbing states the following quantities are character- 
ized by critical exponents | |20| : the total number of 
chaotic sites N{t), the survival probability P{t), the 
mean-squared deviation B?{t) of the turbulent activity 
from the "seed" and the density n{t) of chaotic sites 
within the spreading patch of turbulence. It is assumed 
that they behave according to 



N{t) - 1"- P{t) ~ R^{t) ~ n{t) ~ r 



(7) 



For probabilistic systems it has been conjectured and ver- 
ified numerically pT| , p2| , that the dynamical exponents 
satisfy the generalized hyper-scaling relation rjs-\-5 + 9s = 
dzs/2 where d is the spatial dimension. For systems with 
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FIG. 7: The average spreading of active seeds in the "natural-initial-state" close to criticality. The densities are obtained by 
averaging over 10* realizations, and for clarity we have included three snapshots of this average activity in the bottom rows of 
this figure. Parameters are r — 3 (for all runs), and b = —0.1 (a,d), b — 0.0 (b,e) and b = 0.2 (c,f). The respective values of 
natural-initial-state are 1.170, 1.212 and 1.235. 



a single absorbing state, including DP, one finds that 
S — 9s = and Zg ~ 2/z, reducing the hyper-scahng 
relation to 4(5 -I- 2rjs — dzg. 

Systems with infinite numbers of absorbing states have 
been studied carefully recently and it has been found that 
they differ from the classical ones (with a single absorb- 
ing state belonging to the DP universality class) precisely 
in the non- universality of the spreading exponents [ p2[ . 
Only exponents characterizing quantities averaged over 
surviving runs alone are found to be universal. This im- 
plies that Zs, the sum rjs + 5 and 6s are expected to be 
universal, while 77^ and 5 individually are not. Only for 
the so called "natural-initial-state" are the DP values 
found for the exponents characterizing quantities aver- 
aged over all runs. Such a particular state is constructed 
by letting the system evolve at criticality from homoge- 
neous initial conditions, where all sites initially are in the 
active phase, until an absorbing configuration is reached. 

After a few spreading experiments in our CML we in- 
deed observed that the propagation of activity from the 
initial seed through the laminar region depended strongly 
on the configuration of the laminar state surrounding 
the seed. Moreover the dynamical exponents varied with 
this configuration, thus being non-universal. So the non- 
unique absorbing state of our CML (any configuration 
with all M- values above unity and w-values not too large 
will be absorbing) leads to behavior as can be expected 
for DP with an infinite numbers of absorbing states. 
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FIG. 8: Average spreading activities near criticality for r = 
3, fe = 0.2 (a) and r = 3, 6 = (b). In comparison to the 
spreading into the "natural-initial state" as shown in Fig. ^, 
the value of u in the laminar background has been lowered 
from 1.235 to 1.225 in (a), thus strongly enhancing solitons, 
and has been increased from 1.212 to 1.22 in (b), thus strongly 
suppressing soliton activity. 



We determined the "natural-initial-state" by iterating 
systems of up to 4096 sites from homogeneous initial con- 
ditions until an absorbing configuration is reached. The 
average value of all sites is then used as the value of the 
laminar background. For r — 3.0 we have calculated 
these as 1.235 for b = 0.2, 1.212 for 6 = 0, 1.170 for 
b = -0.1 and 1.0395 for b = -0.2. 

In Fig. ^ we display the average spreading for r = 3 and 
b = -0.1,0 and 0.2. Clearly, for b = -0.1 and for the 
Chate-Manneville model at r = 3, 6 = it is basically 
impossible to estimate the spreading exponents at the 
"natural-initial-state" , since the spreading is dominated 
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0.16495(2) 
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0.16205(1) 
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0.17(1) 
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0.16368(2) 
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0.20(1) 




-0.1 


0.35203(1) 


1.52(3) 


0.02(2) 




0.0 


0.35984(3) 


1.42(2) 


0.18(1) 




0.1 


0.3393(1) 


1.48(2) 


0.155(1) 




0.125 


0.34745(5) 


1.53(2) 


0.15(1) 




0.15 


0.35680(5) 


1.57(1) 


0.159(3) 




0.175 


0.36545(1) 


1.58(1) 


0.16(1) 




0.2 


0.37323(1) 


1.58(1) 


0.16(1) 


DP 






1.58074 


0.15947 



TABLE L The critical exponents z and 6 — l3/v\\ for our 
CML. The values for DP are taken from (l|]. 



by the solitons (see Fig. |7b,e). This behavior is distinctly 
different from what is observed in the various systems 
belonging to the DP universality class. We have therefore 
only estimated the spreading exponents for b — 0.2 at 
branch A, and b = —0.2 and branch B; in both cases the 
solitons are not dominant. 

Note that the strength of the spreading solitons can 
be altered by changing the value of the laminar back- 
ground. By increasing the background to a value above 
the "natural-initial-state" one, the solitons can be sup- 
pressed, while they can be enhanced by a decrease of the 
laminar value; see Fig. |[ 

Our estimates of the dynamical exponents have been 
done for simulations with a maximum time of 2000 it- 
erations. An active seed is placed in the center of the 
lattice, surrounded by a laminar background. The seed 
consists of two active sites, each of which is assigned a 
random number in the chaotic regime, such that the lo- 
cation is fixed but the values of the active sites differ for 
each trial in the ensemble. The ensemble size Ns used 
for statistical averaging and the number of sites in the 
lattice L have been adjusted to the number of surviving 
runs the different setups produced, and how far the tur- 
bulence propagated out from the seed. A minimum of 
200 surviving runs have been used in the averaging. 

Our results in table || agree rather well with previously 
obtained results for probabilistic systems with an infinite 
number of absorbing states. In particular, the exponents 
averaged over surviving runs alone definitely seem to be 
universal as long as the background does not deviate to 
much from the "natural-initial-state" one. While the val- 
ues for the sum S + rjs are very close to the DP value of 
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1.26523 


0.15947 


0.31368 


0.15947 





TABLE IL Estimated spreading exponents for r — 3.0 for 
background values Xi. The deviation from the hyper-scaling 
relation for d = 1 is defined as A = Zs/2 — {rjs+S + 6s). Note 
that for b = 0.2, the natural background state has Xi « 1.235; 
for this value, the exponents 5 and rjs are close to their DP 
values. Similarly for b — —0.2, the natural background state 
has Xi ~ 1.1395; again S and ijs are close to their DP values 



0.47315(7), our results for Zg and 9s deviate from their 
respective DP values. A very interesting observation is 
that the hyper-scaling relation is satisfied (A ~ 0) for 
the majority of different background values. Only for 
the highest values are significant deviations encountered. 



III. STOCHASTIC MODEL. 

The propagating structures, which are observed in the 
CML that we studied in the previous section, appear to 
play an important role for the transitional behavior. It 
is, however, numerically very demanding to obtain good 
statistics for large CMLs and long times. As pointed out 
already in the introduction, this is the reason why one 
tries to map such deterministic models to simple stochas- 
tic models. Not only may there be more hope to under- 
stand such models analytically, they also are much easier 
to handle from a computational point of view. 

In this section we will introduce and study a very 
simple extension of the Domany-Kinzel cellular automa- 
ton which itself is a simple model showing DP behavior. 
While for the Domany-Kinzel automaton every site can 
only be active or inactive, we will allow sites to either 
contain a left or right traveling soliton. As in the CML, 
these solitons should be generated from active sites only, 
and we wish to be able to tune their typical lifetime. The 
only process in which these solitons aid the spreading of 
activity is by collisions: for simplicity we assume that 
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(i) (li) 
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(III) 
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(IV) 
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FIG. 9; Schematic definition of our stocliastic model. In 
this picture, the bottom two circles denote possible incoming 
states, the circle at top denotes the possible stochastic out- 
come. Our model is defined on a diamond lattice and so one 
only needs to define the probabilities for certain offspring (ac- 
tive, inactive, right or left traveling soliton) as a function of its 
two predecessors. Empty circles depict inactive states, black 
circles are active, "R" and "L" denote right and left-moving 
solitons, "S" denotes a soliton of arbitrary direction, and "X" 
finally represents any state (L,R, inactive or active). The no- 
tation m(o, b) denotes the minimal value of a and b, and P 
denotes the conditional probability that this out-coming state 
occurs. See text for a more elaborate explanation. 



with probability one a pair of colliding solitons yields a 
singe active site. 

Belovif we will first discuss the definition of our model 
in section 



L III A , We will then d iscuss the mean field equa- 
tions for our model in section III B| , and these will show 
a transition from second to first order behavior. We will 
study the statistical properties of our model in section 
III C . We will illustrate the role of solitons in direct sim- 



ulations of this model; these simulations will point to the 
relevance of large "holes" which cannot be "healed" by 
the solitons. We will discuss the statistical properties 
of our model near the transition from inactive to active 
states in the soliton-dominated parameter regime. We 
will find that the transition is no longer in the DP uni- 
versality class, since no asymptotic scaling regime can be 
reached. While the transition shows some characteristics 
of a first order transition (dependence on initial state for 
example), the asymptotic situation is not entirely clear: 
rather we find a regime of long lived transient states be- 
tween active and inactive regimes. 



A. Definition of model 

The (1+1) dimensional Domany-Kinzel cellular au- 
tomaton is defined on a diagonal square lattice where 
each site can either be active or inactive. The model 
evolves by parallel updates according to the transition 
probabilities p and q, corresponding to the probabilities 
that an empty plus an active site or two active sites 
respectively produce a single active site. The choice 
q = p(2 —p) corresponds to a realization of directed bond 
percolation jl]. 

In our extension the active sites behave like usual di- 
rected bond percolation except from the fact that with 
probability c they can emit a left or right-moving soliton. 
These solitons have a tunable lifetime and travel ballisti- 



cally. We assume that the solitons cannot, by themselves, 
create chaos, except when two solitons collide. 

The updating rules are illustrated in figure ||, where 
the sites can be either inactive (empty), active (black), 
or contain a left (L) or a right (R) moving soliton (S). 

(i) The inactive state; two inactive sites always yield an 
inactive site. This property ensures that there is a unique 
absorbing state. 

(ii) Soliton propagation: a right moving soliton (R) ei- 
ther dies with probability d, or propagates with proba- 
bility (1 — d) when the "O' state to its right is inactive 
or another right-moving soliton. The rule for left-moving 
solitons follows by left-right symmetry. 

(iii) Soliton collision: when two oppositely propagating 
solitons collide, they generate an active site with proba- 
bility one. This is the only process where solitons lead 
to spread of active sites. In principle we could generate 
active sites with a probability less than one, but it may 
be expected that this does not change the behavior of the 
model in a qualitative sense. 

(iv) Single active sites: a single active site, where X can 
either be a soliton or inactive site, leads with probability 
p to a new active site. Note that the spreading of activity 
is thus not enhanced by individual solitons 

(v) Transformation: a single active site can give rise to a 
soliton (S) with probability min{l —p,c); c denotes the 
creation rate of solitons. Such a new soliton can be either 
left or right-moving with equal probability. 

(vi) Pair of active sites: two active sites create a new 
particle with probability q; we restrict ourselves to bond- 
directed percolation and take q — p{2 — p). 

(vii) Soliton creation from pair of active states: similar 
to case (v), a pair of active sites can give rise to solitons 
with probability min(l — q,c). 



B. Mean Field Equations. 

To interpret the physical properties of our cellular au- 
tomaton, a crude insight can be obtained by applying 
Mean-Field Theory. In this approximation it is reason- 
able to ignore the differences between left and right trav- 
eling solitons, and so our mean field equations are for two 
concentrations, those of chaotic sites c and solitons s. 

b. Equation for Chaotic Sites: Chaotic sites can 
emit solitons and can be generated by collisions of two 
solitons; apart from these two rules they behave like DP. 
Thus, without the solitons the rate equation (without 
noise) would be c = bic — b^c^ To incorporate the 
creation of an active site when two solitons collide accord- 
ing to rule (iii), the term &2S^ needs to be added to this 
equation. There is no source term linear in s in the rate 
equation for c, reflecting that we assume that individual 
solitons do not give rise to activity. 

c. Equation for Solitons : There are four processes 
that influence the solitons. Solitons may decay sponta- 
neously according to rule (ii), and this yields a term —a^s 
in the rate equation for s. Solitons also die upon collision 
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FIG. 10: Dynamical system analysis of the mean field equations (|lO[)-(|ll[). The full and dashed curves show the location in 
r, z space of the saddle-node and transcritical bifurcations respectively. The five insets (a-e) show schematically the flow in the 
various regimes of the mean field equations. For more details see text. 



leading to a term cx — s^. Depending on the lifetime of 
the solitons, either of these two terms may dominate and 
so we keep both of them; we will see below that this will 
indeed be a crucial ingredient. Solitons are created from 
active sites according to rule (v) and (vii). While this 
in principle yields source terms in the rate equation of 
s proportional to both c and (? , we only keep the linear 
term, since the prefactor for both these terms will be of 
the same order. Inclusion of the quadratic term does not 
affect the qualitative dynamics. 

The rate equation for the solitons and chaotic sites can 
then be written as 

s = a\c — a2S^ — a-iS (8) 
c = bic + b2s'^ - bac^ (9) 

where the lifetime of the solitons is set by l/as and the 
spreading rate of the chaotic patches by bi. 

These two equations can be simplified by the intro- 
duction of a rescaled of time r and densities S and C to 



be 

S = C-S'^-aS (10) 

C = tqC + S^-uC^ (11) 

where 

"203 anb, ^2&3 

a= — — , ro = ^4^-, u= — 5- (12) 
t^^r, s = ^S, c^^C (13) 

ai02 °2 0,2 

We will now analyze the possible transitions in the 
mean field Eqs. ([l^p^. 

d. Fixed points The fixed points {S*,C*) of the 
rescaled equations (|l^)-(0) satisfy C* = S*^ + aS*, 
where S* is given by solutions to the fixed point equation 

Sf{S) = (14) 
f{S) = uS^ + 2uaS^ + {ua^ -l-r)S- aro (15) 

Apart from the trivial fixed point {S*,C*) = (0, 0) there 
may be either 1 or 3 other fixed points which can be found 
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from solving Eq. (|l5|). It can be shown that Eq. ( |l5| ) 
always has one solution for large negative S. This fixed 
point can be ignored since only points where both S* 
and C* are positive are relevant for our mean field equa- 
tions (remember that S and C are both concentrations). 
The two non-trivial fixed points (5*1 , Cj*) and (52, C2 ) are 
born in a Saddle-Node bifurcation when the discriminant 
of Eq. ( |l5|) becomes negative. Introducing the parame- 
ter z := aru and performing the tedious standard algebra 
yields that this occurs when 

z2-(2-5ro-r2/4) + (l + ro)3=0 (16) 

and so the locus of the saddle-node bifurcation only de- 
pends on rg and z. It can also be shown that at = 
always one of the non trivial fixed points crosses through 
the fixed point of the origin in a transcritical bifurcation. 
The various types of flows that occur as function of z and 
To are illustrated in Fig. ^ 

As shown in Fig. |l^, there are essentially four quali- 
tatively different types of flow and two bifurcations oc- 
curring. We will here discuss these flowtypes and their 
relevance for the dynamics, (a) Only the trivial fixed 
point is present, and is stable. Hence all initial condi- 
tions flow to the absorbing state, (b) For small soliton 
lifetime z > 1, the two non-trivial fixed points {Sl,C'i) 
and (82,02) that are born in a saddle- node bifurcation 
do not lie in the first quadrant and are therefore not rel- 
evant for the mean field equations. Hence the situation 
in (b) means that there is a single relevant fixed point 
at the origin and so the system is in the absorbing state. 
(c) When, for z > 1, tq crosses through zero from be- 
low, {SI, C*) crosses through the origin in a transcritical 
bifurcation. All initial conditions in the first quadrant 
fiow now to {Sl,Ci); the mean field equations indicate 
that there is a finite activity, whose value grows approx- 
imately linearly in tq. The transition at ro = corre- 
sponds to the standard DP transition for z > 1. (d) For 
long soliton lifetimes (z < 1) the two non-trivial fixed 
points {Si,C*) (square) and (5*2, C2) (triangle) are also 
created in a saddle-node bifurcation; but in contrast to 
case (b) both lie in the first quadrant and are therefore 
relevant for the dynamics. Depending on initial condi- 
tions, the final state can either be absorbing or active; 
the incoming manifold of the saddle point acts as a sepa- 
ratrix. The transition that occurs here as the saddle-node 
bifurcation is crossed leads to a finite jump in the value 
of c in the active state, which is indicative of a first or- 
der transition, (e) When, for z < 1, tq crosses through 
zero from below, (S2 , C'2 ) crosses through the origin in a 
transcritical bifurcation. All initial conditions in the first 
quadrant flow now to (<5'2 , C2). 

To study the phase transition we shall primarily vary 
ro while keeping a and u fixed. There are three generic 
choices for z relevant here: 

• z — > 00. In this case the solitons have probability 1 
to die once they are generated, and so the system 
is effectively soliton-free. This is the case of pure 



DP, and the transition takes place at rp = 0. There 
is no hysteresis. 

• z > 1. This is the regime of short soliton lifetimes. 
Here the solitons do not contribute to any change in 
the qualitative behavior. An attractive fixed point 

S = Si emerges for small, positive rp. This 

2 

corresponds to C — ^ , such that this fixed 
point converges towards the DP value CJ' —>■ rg/u 
for large a. As ro — > this fixed point converges 
towards the origin and it changes stability at tq = 
(see figure |l^), implying that the transition is 
continuous. Thus, the transition for small soliton 
lifetimes (z > 1) still takes place at ro = and 
resembles DP. 

• z < 1. This the soliton dominated regime where 
a completely different scenario occurs. For rp > 
the behavior is determined by the stable node at 
5*2 ~ a(z~^/^ — 1). When tq becomes negative this 
fixed point remains stable and away from the origin. 
Simultaneously the origin becomes attractive and a 
saddle appears close to the origin at S — S^ . 
Initial conditions close to the origin will evolve into 
that point, while initial conditions above the stable 
manifold for the saddle located at {S^, C*) will con- 
verge toward the node {S2, C2 ). This will go on un- 
til the saddle {Sl,Ci ) and the node {S2, Cj) merge 
in a saddle- node bifurcation at rg = ^^(z). Below 
this critical point the origin is globally attractive 
and every trajectory in the phase space converges 
towards this. Going back and forth along scenarios 
(a),(d) and (e) there is hysteresis and so for z < 1 
we clearly observe a first order transition. 

For infinite soliton lifetimes (a = 0, z = 0) the critical 
point is shifted down to tq = —1. Setting a = into 
equation (^5|) yields the fixed point S* = ■\/(l + ro)/u 
which shows that the transition is continuous, but with 
(3=1/2 instead of the DP mean-field value, Pdp = 1- 

Finally, at the tri-critical point (z = 1) equation ( p^ ) 
is reduced to a'^f{S) = S^ + 2aS^ - a'^rgS - a^rg. At 
ro = the only non- negative root is 5 = 0, but for small 
positive To a new root appears at S* « a^/ ro/2. The 
transition thus remains at rp = and is continuous, but 
again with /3 = 1/2 instead of the DP value (3dp — 1. 

C. Phenomenology and statistical properties of the 
stochastic model 

Let us now discuss the properties of the full stochastic 
model based on direct numerical simulations. For small 
but finite values of the soliton lifetime {d ^ 0) or for 
sufficiently small production of solitons (c ^ 1) the tran- 
sition from inactive to active states that occurs when p 
is increased is of second order and indeed appears to be 
in the DP universality class. There is, however, also a 
regime in which the model, at first glance, appears to 
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FIG. 11: (a) Large scale dynamics of our stochastic model for d = 0.01, c — 0.1 and p — 0.614. The grey scale corresponds to 
the number of solitons and active sites coarse grained in a cell of 20 space and 20 time units, (b) Blow up of the dynamics 
shown in the top left corner of (a). 



display a first order transition. In the remainder of the 
discussion on the stochastic model we will focus on this 
regime, which shows some interesting new features. 

The phenomenology of this regime will be illustrated 
following Figs. |ll-12 where different aspects of the dy- 
namics of our model are shown. The parameters cho- 
sen are somewhere in the transitory regime, which in the 
mean field description corresponds to the regime with 
two stable fixed points (Fig. [lO|b). 

In Fig. pT| a we show the evolution of our model, start- 
ing from a fully active state. Fig. |ll]b is a close-up of 
the top left corner of Fig. pT| a which shows the dynamics 
of active sites and solitons in detail. At first glance the 
clusters of activity look extremely similar to the ordinary 
Domany-Kinzel Cellular Automaton, but after closer in- 
spection it becomes clear that colliding solitons generate 
new active clusters (one example can be seen in Fig. pl] b 
for n w 50,i « 4825). 

While Fig. |l^ shows both solitons and active sites, we 
have shown the coarse grained activity and solitons sepa- 
rately in Figs. Oa-b. Clearly, the soliton density is more 
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FIG. 12: Concentration of active sites (a) and solitons (b), 
for the state shown in Fig. |l^. 



uniformly spread, and one can think of the coarse grained 
dynamics as active clusters surrounded by clouds of soli- 
tons. 

To gain insight in the statistical properties of our 
model, we have studied the decay of the number of active 
sites as a function of time, for a range of system sizes L 
and parameter values p. Here and in the remainder, we 
keep the soliton-parameters d and c at values 0.01 and 
0.1 respectively. In Fig. |l^ we show the results of these 
calculations for p ranging from 0.612 to 0.621. 

Clearly, the decay of activity looks quite different from 
DP. If we focus on the activity as a function of p for a 
fixed large time r, we find a very abrupt transition at 
p = Pi from an inactive to an active state, with a value 
of the activity given by the "plateau" that can be seen 
in Fig. |l^. This behavior is indicative of a first order 
transition, consistent with the mean field theory for large 
soliton lifetimes. 

Let us focus on the decay of activity in more detail. 
When p is small enough, the activity will decay faster 
than a power-law. When p large enough so that active 
states will spread, we are certainly above the transition 
and the activity will reach a plateau value. For DP, there 
is a critical value where the activity decays as a power- 
law, but this is not necessarily the case in our model. As 
shown in Fig. ^ there is a transient decay towards the 
plateau value which can look like a power-law, as shown 
in more detail in Fig. |l^. For a transient period that in 
this case goes up to t w 10'^, it is possible to find values 
of p such that the decay of m appears to be a power- 
law with a non-DP exponent. For the example shown, a 
reasonable scaling can be obtained over 2 decades. We 
speculate that this may be the origin of the non-universal 
power-laws observed in coupled map lattices Q . How- 
ever, for this to be real asymptotic scaling, one should be 
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FIG. 13: Decay of average activity m for c = 0.1, d = 0.01, 
and p = 0.610, 0.611, . . . , 0.621 (increasing p leads to an in- 
crease of activity; the curves with p = 0.610 and 0.620 are 
thicker). Averages are taken over (a) 2000 systems of L=200, 
(b) 200 systems of L=2000 (c) 20 systems of L=20000. 



able to have this scaUng extend to arbitrary large times; 
however the activity curves for sizes 200, 2000 and 20000 
all bend downwards at nearly the same time; hence there 
is no hope that increasing the systemsize extends the time 
interval over which apparent scaling can be found. For 
long times the activity either decays rapidly, or first hits 
a plateau. Clearly, the transition is not an ordinary 2nd 
order transition. 

Is this transition now an ordinary 1st order transition? 
When we increase the time scale r where we inspect the 
average activity, the transition becomes sharper, and the 
transition value pi will shift. We do not see evidence for 
Pi going to a well-defined asymptotic value when r goes 
to infinity, at least not on the timescales that we can 
probe numerically, and in this sense the transition is not 
truly first order. 

We will argue now that the plateau indeed does not 
represent the truly asymptotic behavior. Let us return 
to Fig. where the activity appears to arrive at the 
plateau (the overall activity appears to approach a con- 
stant). However, around n « 4000, z « 3000 a large 
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FIG. 14: Average activity m for c = 0.1, d — 0.01 in an 
ensemble of 20 systems of L = 20000, showing the appearance 
of quasi-power-law decay. The three curves correspond to 
p = 0.612, 0.616, and p = 0.62 respectively, and are shifted by 
30% for clarity. The three straight lines corresponds to power- 
laws with exponents —0.19, —0.21 and —0.23. In particular 
the scaling in the system for p — 0.612 looks rather convincing 
with an exponent of —0.23. 



"hole" opens up. Once the size of this hole becomes 
larger than twice the lifetime of the solitons, it becomes 
unlikely that colliding solitons will create new activity 
there and "heal " the hole. In fact, for this particular 
example the hole did spread out and the system decayed 
to the inactive state. A closer inspection of the dynam- 
ical states that occur when the activity drops below the 
plateau value, shows that this is the general scenario: the 
nucleation and subsequent spreading of a large inactive 
droplet is what dominates the asymptotic decay of the 
active states here. 

The spreading or shrinking of a large hole is governed 
by the propagation velocity of the coarse grained domain- 
walls between active and inactive patches. What we have 
found is that the apparent first order transition point pi 
does not coincide with the point p* where this velocity 
changes sign. We illustrate this property by following the 
dynamics of a large inactive droplet for p — 0.618, where 




1 n 5000 1 n 1000 



FIG. 15: Dynamics for c = 0.1, d = 0.01 and p = 0.618, 
showing that a hole of size 500 is healed (a) while a larger 
hole of size 2500 is not healed (b) . 
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10 100 1000 i 10000 

FIG. 16: Evolution of the average activity m in 20 systems 
of size 20000 for c = 0.1, d = 0.01 and p = 0.621, i.e., well 
in plateau regime. Here we vary the initial concentration of 
active states by randomly distributing active sites through 
our lattice for i = 1. These activities are, respectively, 
1.0, 0.85, 0.7, 0.55, 0.4, 0.25, 0.15, 0.1, 0.05 and 0.025. For the 
latter two cases, the plateau is not reached, even though ini- 
tially the activity is increasing. That the long time behavior 
(i > 200) depends on the initial concentration is reminiscent 
of a first order transition. 



the system has a well defined plateau in the activity (see 
Fig. |l^). Thus p > pi, but nevertheless, a hole of size 
5000 grows as can be seen in Fig. |l5|. A hole of size 500 
is healed for these same parameter values, as shown in 
Fig. [l|. 

The solitons in our model thus introduce an additional 
length-scale of order 1 /d, and as a result the size of the in- 
active holes becomes of importance. When such holes are 
smaller than 2/d, solitons will penetrate deeply enough 
to collide and create new patches of activity within this 
hole. Since the first order transition is triggered by the 
solitons, this indicates that the velocity of the domain 
wall is not strongly dependent on these degrees of free- 
dom. This is in accord to the results from the mean field 
analysis, if we assume that the properties of the domain 
walls are still governed by the trivial fixed point with 

The difference between the spreading of a small ac- 
tive cluster and the behavior of an homogeneously ac- 
tive state indicates that the initial concentration of active 
sites plays a role. This is illustrated in Fig. ^ where we 
follow the evolution of the activity for a range of initial 
concentrations of activity for p = 0.621. For initial activ- 
ities in the range from 1 to 0.1, the same plateau value 
is reached, but for initial activities of 0.05 and smaller, 
there is an initial increase of the activity after which the 
activity rapidly decays; the plateau is never reached. 

In Fig. |l^ we show the evolution of the activity m di- 
vided by the number of surviving clusters for the same 
parameter values. For small systems these plots are 
very different from the ones averaged over all systems 
(Fig. |l^; in the present case there is a typical activity 
in each system which rapidly disappears. We interpret 
this as further evidence that the nucleation of large holes 
dominates the eventual decay. For larger systems this 




FIG. 17: Average activity divided by the number of active 
systems, for the same parameter values as shown in Fig. |l^. 

effect disappears because the time it takes for a hole to 
engulf the whole system is large. 

IV. DISCUSSION 

The overall picture that emerges from our study is 
that the transition to spatio-temporal intermittency is 
strongly influenced by coherent ballistically traveling 
"solitons" , which, even though they have a finite life- 
time, change the nature of the transition and can intro- 
duce first order like behavior. That such a scenario is 
relevant is supported by recent evidence for a discontin- 
uous transition to spatio-temporal chaos in the damped 
Kuranioto-Sivashinsky equation , which is well-known 
to support localized ballistically moving excitations, or 
"pulses" ||. 

We build our conclusions upon an extension, using two 
dimensional local maps, of the Chate-Manneville coupled 
map lattice. We thereby gain an additional parameter, 
which turns out to tune the importance and life-time 
of the solitons. For this coupled map lattice we find, 
depending on parameters, evidence for both continuous 
phase transitions in the universality class of Directed Per- 
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eolation with infinitely many absorbing states and for 
first order behavior. 

To understand this behavior, we have developed a 
stochastic model generalizing the Domany-Kinzel cellu- 
lar automaton. In this model, the active sites can emit 
solitons and by colliding, the solitons can create new ac- 
tive sites. Simulations of this model, together with the 
appropriate mean field theory, support the existence of 
both continuous and discontinuous transitions. With the 
stochastic model one can look at the behavior on much 
larger length and time scales. One thereby discovers that 
the active states close above the discontinuous transitions 
are actually metastable, and will finally decay when a suf- 
ficiently large droplet nucleates. Such breakdown of first 
order behavior has been predicted by Hinrichsen . We 



have not been able to pin down exactly how the nucle- 
ation time and eventual decay depends on system size, 
but it is not ruled out that the truly asymptotic behavior 
is still given by Directed Percolation in accordance with 
the predictions of Hinrichsen. 

Even though the first order behavior is probably not 
truly asymptotic, it appears very clearly over a surpris- 
ingly long range of intermediate time scales, and would 
thus be relevant in the interpretation of experiments. We 
further show that this feature can lead to long power like 
transients displaying non- universal "critical exponents," 
and we believe that such transients are the origin of 
the observed non-universality in the transition to spatio- 
temporal intermittency. 
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